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Ab initio quantum chemical cluster calculations of Heisenberg exchange in 
linear cuprate antiferromagnets are reported. For Ca2Cu03 I find J a =-136 
meV for the intrachain coupling and J b =-0.52 meV for the interchain 
coupling. For Sr 2 Cu03 I estimate J a =-136 meV and calculate J b =-0.28 
meV. With these values, expressions resulting from a mean field treatment 
of the interchain coupling reproduce the experimental values for the 
staggered magnetic moment and the Neel temperature T N surprisingly 
well. The linear dependence of T N on J b suggests that Neel order persists 
for any finite value of J b . 
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The one-dimensional (ID) spin 1/2 Heisenberg antiferromagnet is a nontrivial, but 
nonetheless integrable [1,2] quantum mechanical model with vanishing long-range 
magnetic order in the ground state. In real systems that contain antiferromagnetic chains 
there always exists interchain coupling. The linear cuprates Ca 2 CuC>3 and Sr 2 Cu0 3 [3] 
presently attract much attention because they have exceptionally small ratio's of inter- to 
intrachain coupling, J b /J a , and thus are nearly ideal ID spin 1/2 systems. This was 
inferred from the magnetic susceptibility and muon spin relaxation (|lSR) experiments 
[4,5]. These materials therefore provide a unique opportunity to compare theoretical 
predictions with experiment, provided that the values for the coupling parameters are 
accurately known. Unfortunately, the values for J a as derived from experiments vary 
considerably. From an analysis [6,7] of magnetic susceptibility measurements [8] a value 
of -J a ~135-160 meV was found. A second susceptibility study obtained a value of 
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J a =-190 meV [9]. An even larger value of J=-260 meV was estimated from the analysis 
of the mid-gap infrared absorption spectrum [10]. For comparison, in the planar cuprates 
J a varies between -100 and -140 meV [11]. For the interchain coupling J b no 
experimental value is available, but mean field expressions of the staggered magnetic 
moment [12,13] and the Neel temperature [12] in terms of J b and J a have been proposed. 
Recent (iSR and neutron scattering experiments [5] give precise values for the staggered 
moment and the Neel temperature of Ca 2 Cu03 and Sr 2 CuC>3. Clearly, accurate theoretical 
estimates for J a and J b would resolve the controversy over J a and permit a check on the 
validity of the mean field approximation of the interchain coupling. 

In this Letter values for J a and J b are computed in a manner that is entirely independent 
from experiment. The computational approach that can achieve this reproduces the in- 
plane exchange coupling in quasi-2D cuprates within experimental accuracy [11,14]. Its 
application to the linear cuprates and extension to the description of interchain coupling is 
straightforward. The approach is an accurately balanced calculation of singlet-triplet 
splitting of a cluster containing two Cu 2+ ions with quantum chemical all-electron 
methods. Using the computed values of J a and J b a detailed comparison of the various 
theoretical predictions and experimental results is made. 

The structure of Ca(Sr)2Cu03 consists of two interpenetrating sublattices of infinite 
chains of S=l/2 [CUO3]" 4 units separated by Ca(Sr) +2 ions [8,15], that are magnetically 
decoupled by symmetry. The magnetic behavior of a single sublattice can be described by 
the spin hamiltonian 

—> — > -> -> 

H = -%c(Ja S ijk* S i+ljk + J b S ijk* S ij+lk)> C 1 ) 

where the summation runs over lattice vectors (ijk). The chains are oriented along a and 
the Cu03-units lie in the a-c-planes, so that J a and J b are the intra- and interchain 
coupling, respectively. 

The computational approach for the intrachain coupling J a follows that of Refs. 
[11,14]. As a first step the (MC)SCF singlet and triplet ground states of the planar, D 2 h 
symmetric, CU2O7 cluster that occurs in the chains are calculated. It is necessary [11,14] 
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to introduce local exchange and correlation effects on the bridging oxygen atom. This is 
done in a second step of the calculation through admixture of excited (MC)SCF states, 
that differ from the MCSCF ground states by an O -> Cu electron excitation. This is an 
example of non-orthogonal CI (NOCI) [16]. NOCI involves the computation of 
Hamilton and overlap matrix elements between determinants constructed from non- 
orthogonal orbital sets [17]. It has the advantage over conventional CI methods that it 
leads to a short, physically transparent wave function and that it is free of size 
consistency errors. 

For the calculation of the interchain coupling J b , I perform all-electron (MC)SCF 
calculations of the singlet and triplet ground states of D 2 h symmetric Cu 2 8 and 
Cu20gM 16 (M=Ca,Sr) clusters (Fig. 1) consisting of two Cu0 4 units that are translated 
with respect to each other by one unit vector along the b-axis, so that they are located in 
next-nearest neighbor chains. Since there are no ligands, the exchange originates simply 
from the overlap between the Cu d-orbitals and J b can be found directly from the 
(MC)SCF ground state singlet-triplet splitting. For a study of the effect of the 
counterions, calculations were also performed on Cu20 8 M 16 (M=Ca,Sr). 

The clusters are embedded in a point charge environment that accurately represents the 
Coulomb potential inside the cluster region of a surrounding infinite lattice of formal 
ionic charges. The use of the simple point charges at nearest neighbor positions may lead 
to spurious occupation of diffuse orbitals [18], because these orbitals feel the strong 
attractive potential of the Cu 2+ and Sr 2+ point charges, without being repelled by the ion 
core [13, 16]. Therefore the potential due to the ions at the nearest neighbour positions to 
the cluster, is modified to be flat inside a small sphere with an ionic radius [19] around 
these point charges. These modified potentials improve the stability against variation of 
the most diffuse components of the basisset, whereas in the cluster region the potential 
due to the modified charges is identical to that of the bare point charges. The complete 
specification of the background potential is available from the authors on request. 

For Cu and O the same contracted Gaussian basis sets are used as in ref. [11], namely 
a 14sllp6d -> 8s6p3d basis set for Cu and a 9s6p -> 3s3p one for O. It was found that 
J b is sensitive to extension of the Cu basis with diffuse d-functions, as will be discussed 
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below in detail. For Ca I employed a 14s8p -> 5s2p basis set of double zeta quality [20] 
and a 12s6p -> 4s2p minimal basis set [21] and for Sr a (15s,9p,3d) -> (5s,3p,ld) 
minimal basis set [21]. 

First consider the intrachain coupling. Following Ref. [11], in an SCF wavefunction 
for the lowest triplet state of the [Cu 2 7 ]~ 10 model clusters two open shell orbitals occur, 
denoted by d g and d u and transforming as a g and b 3u , respectively. These orbitals can be 
considered even (g) and odd (u) linear combinations of Wannier-like orbitals, dj and d 2 , 
localized at the Cu. The triplet wavefunction can be written as 

Y t = I ood g d u I = I oodid 2 1 (2) 

Here a denotes the 0(2p a ) orbital at the bridging oxygen, which has the same symmetry 
as d u . For clarity all other closed shell orbitals are suppressed in the notation. The singlet 
corresponding to Eq.( 2) can be written as 

= (2+2S 2 y 1/2 I 00 ( d : d 2 - d : d 2 ) I 

= (2+2S 2 y 1/2 {I (1+S) oodgdg - (1-S) aad u d u I }, (3) 

In Eq. (3) d^ and d 2 have an overlap S = <di I d2 >, which constitutes an additional 
variational parameter. As also found for the planar cuprates, the triplet and singlet SCF 
ground states are very well characterised by Cu +2 (3d 9 ) and 0~ 2 (2p 6 ) and the Cu holes 
have almost pure 3d(x 2 -y 2 ) character in a Mulliken orbital population. The singlet states 
(3) with optimized orbitals have S-0.04 and their energies are 20 to 30 meV below the 
triplets. The wave functions (2) and (3) form an appropriate starting point of a balanced 
calculation of the singlet-triplet splitting [22] . At this level of the calculation the splitting 
corresponds to Anderson superexchange [23], which is due to charge transfer excitations 
of the type di) 1 d2) 1 — > di) 2 d 2 )°+di)°d2) 2 that can only occur for the singlet. 

As a next step I admix to *P t and *P S relaxed charge transfer excitations of the form 
di) 1 o) 2 d2) 1 ^di) 1 o) 1 d2) 2 ±di) 2 a) 1 d 2 ) 1 . Admixture of unrelaxed excitations of this type 
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has no effect, but this is quite different for the relaxed charge transfer excitations. As 
discussed in full detail in Ref. [11], the charge transfer states, of 3 B3 U and *Ag 
symmetry, can be written as 



¥ t * = I d u d u d g o I , 



(4) 



(1+S*) I dgdguiuj I - (1-S*) I d g d g u 2 u 2 I 
V(2+2S* 2 ) 



(5) 



In the excited singlet state (5) again an overlap between open shell orbitals is 
allowed, now between the 0(2p a ) and the Cu(3d u ) orbitals. Moreover, the open shell 

orbitals are allowed to mix because they have the same symmetry. It turns out that the 
orbitals uj and u 2 approximately correspond to bonding and antibonding combinations of 
d u and a. It should be noted that in the wave functions (2-5) the orbitals are separately 
optimized. An energy separation of about 10 eV between ground and excited states is 
found. As for the planar cuprates, S*~0.55 is obtained. This result can be interpreted as 
a tendency towards covalent bond formation between the bridging O and the remaining 
Cu 2+ neighbour in the excited singlet state, which is enhanced by orbital relaxation. 
These effects are absent in *P t *. The CU2O7 wave functions of the linear cuprates 
obtained here closely resemble those of the analogous clusters for the planar cuprates. 

The NOCI calculation involves a 2x2 non-orthogonal diagonalization for the triplet. As 
argued in Ref. [11], the singlet involves a 5x5 non-orthogonal diagonalization, since S, 
S* and the parameter determining uj and u 2 as orthogonal combinations of d u and a have 
to be reoptimized. The results are displayed in Table I. For both Ca 2 Cu03 and Sr 2 Cu03 
a value of J a =l 19 meV is obtained, which is close to value for La 2 CuC>4 reported in ref. 
[Van Oosten, 1996 #205]. When the 16 nearest Ca neighbors are included, as all 
electron atoms, the value for Ca 2 Cu03 increases to 136.5 meV. In view of the Cu 2 C>7 
results, it is reasonable to adopt this value also for Sr 2 Cu03. This value agrees with 
previous analyses [6,7] have given of the susceptibility data [8]. Much larger values, 
such as J a =-190 meV [9] and J a =-260 meV [10] should be excluded. 
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Let us now turn to the calculation of the interchain exchange J b from the singlet-triplet 
splitting of the [C^Og]" clusters. There are no bridging O so that J b equals the singlet- 
triplet splitting at the (MC)SCF level and is obtained directly from wave functions 
analogous to equations to (2) and (3). Note that the orbital d u now transforms as bj u and 
that the a orbital can be suppressed in the notation. Since the Cu-Cu distance of 3.278 A 

o 

(Ca2CuC>3) and 3.494 A (S^CuC^) is rather larger than a Cu-Cu inter-atomic distance, I 
investigated the effect of augmenting the basis sets with diffuse Cu d- and O p-functions. 
The effect of extending the d-basis is substantial and saturates after addition of a=0.05, 
0.016 and 0.005 functions. An extra a=0.04 O p-function has negligible effect. The 
results are listed in Table II. The diffuse functions describe the exponential decay of the 
d-orbital in the region of space where the overlap is generated and contribute little to the 
total energy. The interchain coupling saturates towards values of J b =0.33 meV 
(Ca 2 Cu03) and J b =0.23 meV (Sr 2 CuC>3). The overlaps of the d-orbitals are close to 
4 %o. The calculation was repeated for Cu20gM 16 (M=Ca,Sr) clusters that include the 
nearest 16 counterion neighbors, using the extended Cu basis. As seen from Table III, 
the effect of the counterions is a 50 % increase of -J b for Ca2Cu03 and 20 % increase for 
Sr 2 Cu03, where it makes no difference wether the double zeta or the minimal basis set is 
used for Ca. These values can be inserted into expressions for the staggered moment and 
the Neel temperature obtained from a mean field treatment of J b [12,13]. It should be 
noted that a chain couples antiferromagnetically to two out of four next-nearest 
neighbors, whereas the coupling to the other two (J c ) is ferromagnetic and should be 
considerably smaller than J b . The coupling to the four nearest neighbors cancels in the 
ordered state. As can be seen from Table IV, this yields values in excellent agreement 
with experiment [5]. Equally good agreement exists between the measured values and the 
mean field prediction for T N [12], when the calculated J a and J b are used. The situation is 
graphically displayed in Fig. 2. In view of the small number of neighbor chains, it is 
quite surprising that mean field theory should perform so well in the present case. The 
linear dependence T N as a function of J b at J b /J as small as 2.4x10 supports the idea 
that long range magnetic order persists at arbitrarily small coupling ratio [24]. 
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In summary, the intra- and interchain magnetic coupling, J a and J b , in Ca 2 Cu03 and 
Sr 2 CuC>3 were computed . I obtain values of J a =-1 19 meV for the intrachain coupling in 
both Ca2Cu03 and S^CuC^, respectively, with CU2O7 clusters. This value is close to 
the value previously obtained for La 2 CuC>4. For the Cu 2 07M 16 (M=Ca,Sr) clusters J a is 
enhanced to -136 meV, which agrees with several analyses of the magnetic susceptibility 
of Sr 2 Cu03. For the interchain coupling I obtain J b =-0.52 meV and J b =-0.28 meV, 
respectively. With these values mean field expressions predicts magnetic moments and 
Neel temperatures in close agreement with experimental data. This constitutes strong 
support to the mean field treatment of the interchain coupling and to the idea that magnetic 
order occurs for any finite interchain coupling. 

I am grateful to Stephan Drechsler, Walter Stephan and Karlo Penc for stimulating 
discussions. 
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Tables 



a b c 

Sr 2 Cu0 3 (Cu 2 7 ) -32.8 -118.7 -100..-260 1 
Ca 2 Cu0 3 (Cu 2 7 ) -36.5 -118.8 — 
Ca 2 CuQ 3 (Cu 2 7 Ca 16 ) -48.1 -136.5 



Table I. Calculated intrachain coupling J a of the linear cuprates Sr 2 Cu0 3 and Ca 2 Cu0 3 : 
a) Triplet- singlet splitting of Eqs. (2,3); b) Non-orthogonal CI involving (2,3) and (4,5); 
d) Experiment, ^efs. [6-10]. 



added diffuse exponent 
2p Cu 3d 


Jb(meV) 
Sr 2 Cu0 3 Ca 2 Cu0 3 




-0.167 


-0.209 


0.05 


-0.176 


-0.252 


0.05, 0.016 


-0.226 


-0.321 


0.05, 0.016, 0.005 


-0.228 


-0.330 


0.04 


-0.166 


-0.213 



Table II. Calculated interchain coupling J b of Sr 2 Cu0 3 and Ca 2 Cu0 3 from Cu 2 Og 
clusters, for various basis sets. 
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cluster 


Jb(meV) 




2p 


Cu 2 8 


Cu 2 8 M 16 


Sr 2 Cu0 3 


-0.228 




-0.276 


Ca 2 Cu0 3 


-0.330 




-0.522 



Table III. Calculated interchain coupling J b of Sr 2 Cu0 3 and Ca 2 CuC>3 for Cu 2 Og and 
Cu 2 8 M 16 (M=Sr,Ca) clusters. 



magnetic moment (|X B ) Neel temperature (K) 





Mean field 1 ' 2 


(iSR 3 


Mean field 2 


(iSR 










n-diffraction 


Sr 2 Cu0 3 


0.065 


0.06(1) 


5.58 


5.41(1) 


Ca 2 Cu0 3 


0.089 


0.09(1) 


10.12 


10 



Table IV. Mean field and experimental staggered magnetic moment and Neel temperature 
in Sr 2 Cu0 3 and Ca 2 Cu0 3 , using Cu 2 8 M 16 (M=Sr,Ca) clusters. : Ref. [13]. 2 Ref. [12]. 
3 Ref. [5]. 
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Figure captions 



FIG. 1. The Cu 2 8 M 16 (M=Sr,Ca) cluster. 

FIG. 2. Staggered magnetic moment versus interchain coupling IJ b l. The solid curve 
represents the mean field expression. The points correspond to Table IV. 
FIG. 3. Neel temperature versus the interchain coupling IJ b l. The solid curve represents 
the mean field expression. The points correspond to Table IV. 
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